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Abstract. A new solvable many-body problem is identified. It is characterized by nonlinear 

t-H Newtonian equations of motion ("acceleration equal force") featuring one-body and two- 

body velocity-dependent forces "of goldfish type" which determine the motion of an arbitrary 
number A" of unit-mass point-particles in a plane. The N (generally complex) values z n (t) 
at time t of the A" coordinates of these moving particles are given by the N eigenvalues of 

I— j a time-dependent N x N matrix U(t) explicitly known in terms of the 2N initial data z n (0) 

^ and z n (0). This model comes in two different variants, one featuring 3 arbitrary coupling 

constants, the other only 2; for special values of these parameters all solutions are completely 
periodic with the same period independent of the initial data ("isochrony"); for other special 

| | ' values of these parameters this property holds up to corrections vanishing exponentially as 

t — > oo ("asymptotic isochrony"). Other isochronous variants of these models are also 
reported. Alternative formulations, obtained by changing the dependent variables from the 
N zeros of a monic polynomial of degree A" to its N coefficients, are also exhibited. Some 
mathematical findings implied by some of these results - such as Diophantine properties of 
the zeros of certain polynomials - are outlined, but their analysis is postponed to a separate 

y—i paper. 

Key words: nonlinear discrete-time dynamical systems; integrable and solvable maps; 
iy-} isochronous discrete-time dynamical systems; discrete-time dynamical systems of goldfish 
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1 Introduction 



^ The technique used in this paper to identify a new solvable many-body problem has become by 

now standard. Its more convenient version starts from the identification of a solvable matrix 
problem characterized by two first-order, generally autonomous, matrix ODEs defining the time 
evolution of two N x N matrices U = U(t) and V = V(t): 

U = F{U,V), V = G(U,V), (1.1a) 



where the two functions F(U,V), G(U,V) may depend on several scalar parameters but on 
no other matrix besides U and V. Here and hereafter superimposed dots denote of course 
differentiations with respect to the independent variable t ("time"). The solvable character of 
this system amounts to the possibility to obtain the solution of its initial-value problem, 

U(t) = U(t;U ,V ), V(t) = V(t;U ,V ), U = U(0), U = V(0), (1.1b) 



*This paper is a contribution to the Special Issue "Geometrical Methods in Mathematical Physics". The full 
collection is available at http://www.emis.de/journals/SIGMA/GMMP2012.html 
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with the two matrix functions U(t;Uo,Vo), V(t;Uo,Vo) explicitly known. For instance, in this 
paper we shall focus (see (3.3) below) on a system (1.1a) which can be related to a linear 
(third-order, matrix) evolution equation (see (3.3) below), so that its initial- value problem can 
be explicitly solved. 

One then introduces the eigenvalues z n (t) of one of these two matrices, by setting, say 

U(t) = R{t)Z{t)[R{t)}-\ Z{t) = diag[z„(t)]. (1.2) 

Remark 1.1. The diagonalizing matrix R(t) is identified by this formula up to right- multipli- 
cation by an arbitrary diagonal N x N matrix D(t), R(t) R(t) = R(t)D(t). 

Then one introduces a new N x N matrix Y(t) by setting 

V(t) = R(t)Y(t)[R(t)]-\ Y(t) = [R^Vmit), (1.3a) 

where of course the N x N matrix R(t) is that defined above, see (1.2). This matrix Y(t) is of 
course generally nondiagonal: 

(t) + (1 - 5 nm )Y nm (t). (1.3b) 

Notation 1.1. Indices such as re, m, i run from 1 to N (unless otherwise explicitly indicated), 
and N is an arbitrary positive integer (indeed generally N > 2). 5 nm is the Kronecker symbol, 
<5nm = 1 if re = m, 5 nm = if n ^ m. Note that hereafter we use the notation Y nm to denote 
the N (N — 1) off-diagonal elements of the N x N matrix Y. 

It often turns out that the time evolution of the eigenvalues z n (t) is then characterized by 
a system of N second-order ODEs which read as follows: 



N 



f(z n , Z n ) + ^ 



£=lj=n 



v v 9 { ) {z n ,z n )g { ^>{z h z i ) 

Z n ~ Zt 



(1.4) 



where the three functions f(z, i), g^\z, z) and g( 2 \z, z) can be computed from the two matrix 
functions F(U,V), G(U,V) (see below). It is then natural to try and interpret this system 
of ODEs as the Newtonian equations of motion ("acceleration equal force") of a many-body 
problem characterizing the motion of N particles whose coordinates coincide with the ./V eigenva- 
lues z n (t); an iV-body problem generally featuring one-body and two-body velocity-dependent 
forces, with the N(N — 1) quantities Y n (Y~i n playing the role of "coupling constants". But 
these quantities are not time-independent, nor can they be arbitrarily assigned: they are the 
off -diagonal elements of the N x N matrix Y, hence they should be themselves considered as 
dependent variables, the time evolution of which is characterized by the system of N(N — 1) 
ODEs implied for them by (1.1a) via (1.2) and (1.3). 

Two options are then open to provide nonetheless a "physical" interpretation for the equations 
of motion (1.4). 

One option that we do not pursue here is to provide some kind of "physical" interpretation 
for these 7V(iV — 1) quantities Y n gYi n as additional (internal) degrees of freedom of the moving 
particles. 

A second option - the one which we pursue below - is to find a (time-independent) ansatz 
expressing the N(N — 1) quantities Y n £ in terms of the ./V coordinates z m , or possibly of the 2N 
quantities z m , z m ; an ansatz consistent with the N(N — 1) equations of motion satisfied by the 
N(N — 1) quantities Y n £, which satisfies these equations either identically (i.e., independently 
from the time evolution of the N coordinates z m (t)) or, as it were, self- consistently (i.e., thanks 
to the time evolution (1.4) of the N coordinates z m (t) with the N(N — 1) quantities Y n £ assigned 
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according to the ansatz). Given a matrix system of type (1.1a) no technique is known to assess 
a priori whether or not such an ansatz exist. However the experience accumulated over time 
suggests that, if such an ansatz does exist, it has one of the following two forms: 

ansatz 1 : Y n £ = - — - — — — ; (1.5a) 

z n - z e 



ansatz 2 : Y, 



ni = {g {1) (z n )gW(z e ) [z n + fW(z n )] [z t + f< 2 \z £ )] } V2 ; (1.5b) 

of course with the functions appearing in the right-hand side of these formulas chosen appro- 
priately. And as a rule the ansatz 1 should work identically, i.e. independently of the time 
evolution of the coordinates z m (t); while to ascertain the validity of ansatz 2 the equations of 
motion (1.4) should be used self- consistently. 

For instance in the very simple case of the equations of motion (1.1a) with F(U, V) = V and 
G(U, V) = implying U = and U(t) = Uq + Vot, V(t) = Vq, both ansatze exist: the ansatz 1 
reads in this case Y n £ = ig/(z n — Zf) with i the imaginary unit (introduced for convenience) 
and g an arbitrary constant, and it yields the prototypical "CM" model characterized by the 
equations of motion 



N 

2g 2 Y, 0n-^)~ 3 ; (1-6) 

=1,^71 



while the ansatz 2 in this case reads Y n g = (znZp) 1 / 2 , and it yields the prototypical "goldfish" 
model characterized by the equations of motion 

± ,l,a, 

Nomenclature and historical remarks. The model characterized by the Newtonian 
equations of motion (1.6) - which obtain from the Hamiltonian 

N , N 



h C m(z,p) = - J2pI + \a 2 



n=l n,m=l,ny^m 

- is usually associated with the names of those who first demonstrated the possibility to treat 
this many-body problem exactly, respectively in the quantal [1] and in the classical [2] contexts; 
accordingly, we usually call "many-body problems of CM-type" those featuring in the right- 
hand ("forces") side of their Newtonian equations of motion a term such as that appearing in 
the right-hand side of (1.6) (in addition of course to other terms). Such models are typically 
produced by the ansatz 1. 

The solvable character of the many-body problem characterized by the Newtonian equations 
of motion (1.7a) - which is also Hamiltonian, for instance with the Hamiltonian [3, 4, 5] 



N 



H g oidU,p) = Y 



71=1 



N 

exp(p n ) Y\ ( 

m=l,m^n 



- is demonstrated by the following neat Prescription [5, 6, 7]: the iV values of the coordina- 
tes z n (t) providing the solution of the initial-value problem of the equations of motion (1.7a) 
are the iV zeros of the following algebraic equation for the variable z: 

v \ in(0 M = -• (i-7b) 
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(Note that this formula amounts to a polynomial equation of degree N in z, as it is immediately 

N 

seen by multiplying it by the polynomial f\ [ z ~ z m(0)])- In [7] it was suggested that this 

m=l 

model, in view of its neat character, be considered a "goldfish" (meaning, in Russian traditional 
lore, a very remarkable item, endowed with magical properties); accordingly, we usually call 
"many-body problems of goldfish-type" those featuring in the right-hand ( "forces" ) side of their 
Newtonian equations of motion a term such as that appearing in the right-hand side of (1.7a) 
(in addition of course to other terms). Such models are typically produced by the ansatz 2. 

The original idea of the approach described above is due to Olshanetsky and Perelomov, 
who introduced it to solve, in the classical context, the many-body model characterized by the 
Newtonian equations of motion (1.6) [8]. For a more detailed description of their work see 
the review paper [9], the book [10], Section 2.1.3.2 (entitled "The technique of solution of Ol- 
shanetsky and Perelomov") in [11], and other references cited in these books. In Section 4.2.2 
(entitled "Goldf ishing" ) of [5] several many-body problems, mainly "of goldfish type", are re- 
viewed, the solvable character of which has been ascertained by this approach; and several other 
such models are discussed in more recent papers [12, 13, 14, 15]. 

The present paper provides two further additions to the list of solvable many-body problems 
"of goldfish type" ; and we expect that other items will be added to this list in the future, possibly 
by a continuation of the case-by-case search of solvable matrix evolution equations allowing - via 
the route outlined above and described in more detail, in a specific case, below (see Section 3) - 
the identification of working ansdtze leading to new systems of Newtonian equations (thereby 
shown to be themselves solvable, inasmuch as their solution is reduced to the algebraic task of 
computing the N eigenvalues of an explicitly known N x N matrix). The identification of a new 
model of this kind constitutes - in our opinion - an interesting finding (even if several analogous 
models have been previously discovered) : we view these many-body problems as gems embedded 
in the magma of the generic many-body problems which are not amenable to exact treatments 
(although the latter include of course more examples of applicative interest and are also mathe- 
matically interesting to investigate the emergence and phenomenology of chaotic behaviors). 

In the following Section 2 the main findings of this paper are reported, including in particular 
a description of two new solvable many-body problems (one featuring 3, the other only 2, 
a priori arbitrary parameters), of the algebraic solution of their initial- value problems, and of the 
variety of behaviors (including isochrony and asymptotic isochrony) featured by them for certain 
assignments of their parameters; two additional isochronous systems are moreover exhibited in 
Section 2.1. In Section 3 these results are proven. Section 4 provides the alternative formulations 
of these models, obtained by changing the dependent variables from the N zeros of a monic 
polynomial of degree N to its N coefficients. A final Section 5 entitled "Outlook" outlines 
further developments, including in particular the identification of Diophantine properties of the 
zeros of certain polynomials; but their detailed discussion is postponed to a separate paper. 



2 Main findings 



The two models treated in this paper are characterized by the following two sets of Newtonian 
equations of motion "of goldfish type" . 
Model (i): 



z n = -3z n z n + 7i n - z n + jz n + [-a + 6(7 + b)]z n + 0(7 + b) 
^ 2(z n + a + bz n + z n ) (ze + a + bze + zf) 

e=ij^ n L Zn Zl 

where a, b and 7 are 3 a priori arbitrary parameters. 



(2.1a) 
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Model (ii): 

z n = —3z n z n — 3bz n — z\ — 3bz^ — 2 (a + b 2 )z n — 2ab 

2(z n + a + bz n + 4/2) (z e + a + bz e + zj/2) 



N 

+ E 

where a and b are 2 a priori arbitrary parameters. 



(2.1b) 



Notation 2.1. These models describe the motion of an arbitrary number N (generally N > 2) 
of points moving in the complex z-plane (but see below Remark 2.1). Their positions are 
identified by the complex dependent variables z n = z n (t). The independent variable t ("time") 
is instead real. Superimposed dots denote of course time-differentiations. The parameters 
featured by these models are generally arbitrary complex numbers; unless otherwise specified 
when discussing special cases. 

Remark 2.1. It is possible to reformulate these models so that they describe the motion of 
pointlike "physical particles" moving in a real - say, horizontal - plane: see for instance, in [11], 
Section 4.1 entitled "How to obtain by complexification rotation-invariant many-body models 
in the plane from certain many-body models on the line". This task is left to the interested 
reader. But hereafter we feel free to refer to the models identified by the Newtonian equations 
of motion (2.1) as many-body problems characterizing the motion of N particles in a plane. 

Remark 2.2. Additional parameters could be inserted in these models by shifting or rescaling 
the dependent variables z n or the independent variable t. We will not indulge in such trivial 
exercises (see also Remark 3.1 below). 

The solvable character of these two models is demonstrated by the following 

Proposition 2.1. The solution of the initial-value problems of the two many-body models cha- 
racterized by the Newtonian equations of motion "of goldfish type" (2.1) are given by the eigen- 
values of the N x N matrix U(t), the explicit expression of which in terms of the 2N initial 
data z n (0), z n (0) and of the time t is given by the following formulas (see (3.6a) with (3.7)): 

U(t) = i{/ + Aexp[i(uj2 — <x>i)i] + B exp[i(ix>3 — ^i)i]} 1 

x {coil + W2^4 exp[i(w2 — 0Ji)t] + LU3B exjp[i(uj3 — wi)£]}, (2.2a) 

with the two constant N x N matrices A and B defined as follows: 

A = -(wi - w 3 )(w 2 - ^y 1 ^ ~ K w l + UJ 3){Vo + w|) 1 (U - i^y 1 ] 

x [/-i( W2 + W 3)(y + a;|)- 1 ([/ _i a ;3)-i]- 1 , (2.2b) 

B = (u 1 - uj 2 )(uJ2 - [I - + w 2 )(Vo + w|) (U - i^Y 1 } 

x [l-i(uj 2 +u z )(V a + ujly\u -iu 2 y l Y\ (2.2c) 

Here and hereafter I is the N x N unit matrix, the N x N matrix Uq = U(0) is diagonal and is 
given in terms of the initial particle positions z n (0) as follows, 

U = diag[z n (0)], (2.3) 

while the N x N matrix Vq = V(0) is the sum of a diagonal and a dyadic matrix, being given 
componentwise by the following formulas in terms of the initial particle positions z n (0) and 
velocities z n (0): 

(V ) nm = -S nm [a + bz n {0) + (c - 1)4(0)] + V n V m , (2.4a) 
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Vn = [in(0) + a + bz n (0) + cz 2 n (0)] 1/2 , (2.4b) 

with c = 1 for model (i) and c = 1/2 for model (ii) (see (2.6)). As for the 3 constants ujj 
appearing in (2.2), they are defined by the following formula (see (3.5b)): 

uj 3 + ijuj 2 + (3uj — ia = (uj — uj\)(uj — uj2)(uj — UJ3) (2.5a) 

implying 

a = — ioJioj2UJ3, P = oj\U)2 +^2^3 + CJ3W1J 7 = i(oo\ + W2 + W3), (2.5b) 

wit/i a and /3 respectively a, (3 and 7 given by the following expressions for model (i) respectively 
for model (ii) (see (3.21a) respectively (3.21b)): 
model (i): 

c = l, a = a(j + b), p = -0 + 6(7 + 6); (2.6a) 
model (ii): 

c=-, a = -2a6, ^ = -2(a + 6 2 ), 7 = -36. (2.6b) 

Note that the 3 a priori arbitrary parameters Uj have the dimension of an inverse time; above 
and hereafter we assume for simplicity that they are different among each other (except in the 
following Subsection 2.1, where they are all assumed to vanish). 

It is plain (see (2.2)) that, if the 3 constants Uj are integer multiples of a single real constant uj, 

ujj = kjUJ, j = 1,2,3, (2.7a) 
then the matrix U(t) is periodic, 

U(t) = U(t ±T), (2.7b) 
with period 

T = ri- (2-8) 

\uj\ 

Here and throughout the 3 parameters kj are integer numbers (positive, negative or vanishing, 
but different among themselves); their definition, as well as that of the positive parameter uj, is 
made unequivocal (up to permutations; once the 3 parameters a, P, 7 are assigned, compatibly 
via (2.5b) with (2.7a)) by the requirement that this positive parameter uj be assigned the largest 
value for which (2.7a) holds. 

More generally, if the real parts of the 3 constants ujj are integer multiples of a single real 
constant to and the imaginary parts of 2 of them coincide while the imaginary part of the third 
is larger, say, 

Ke(ujj) = kjUJ, j = 1,2,3; Im(tJi) = Im(o;2) < Im(cj3), (2.9a) 

then the matrix U(t) is asymptotically periodic with period T, namely it becomes periodic with 
period T in the remote future up to exponentially vanishing corrections, so that 



lim \U(t) -U(t±T)\ = 0. 

t— >oo 



(2.9b) 



Another New Solvable Many-Body Model of Goldfish Type 



7 



While for generic values of the parameters, implying via (2.5) that the imaginary parts of 
the 3 quantities ujj are different among themselves, 

Im(a;i) ^ Im(cj 2 ), Im^) / Im(o;3), Im^) / ~im.{uj\), (2.10a) 

then clearly U(t) tends to a time-independent matrix as t — > ±oo, 

U(t) -»■ U(±oo). (2.10b) 

t— >±oo 

Let us emphasize that these outcomes, (2.7b) or (2.9b) or (2.10b), obtain - provided the 
3 parameters a, /3, 7 satisfy (2.5b) with (2.7a) or (2.9a) or (2.10a) - for arbitrary initial data, 
and that the period T is as well independent of the initial data. Let us also note that these 
properties hold unless U (t) is singular; clearly a nongeneric circumstance, in which case U (t) 
would in fact still feature the properties indicated above, but only in the sense in which, for 
instance, the function t&n[u(t—to)/2] (with uj and to two real numbers) is periodic with period T. 

These properties of the N x N matrix U(t) carry of course over to its eigenvalues z n (t), 
hence to the generic solutions of the many-body problems "of goldfish type" (2.1); these models 
are therefore isochronous respectively asymptotically isochronous if their parameters satisfy the 
relevant conditions, see (2.5b) with (2.7a) respectively (2.9a). 

Remark 2.3. Let us however recall that the periods of the time evolution of individual eigen- 
values of a periodic matrix may be a (generally small) positive integer multiple of the period of 
the matrix, due to the possibility that different eigenvalues exchange their roles over the time 
evolution (for a discussion of this possibility see [16], where a justification is also provided of 
the statement made above that the relevant positive integer multiple is "generally small"). 

2.1 Two additional isochronous many-body models 

Additional isochronous models obtain by applying, to the special cases of the two many-body 
models (2.1) with all parameters vanishing, 



2(z n + czl){zj + czj) 

Zn — Z£ 



(2.11) 



with c = 1 respectively c = 1/2, the standard "isochronizing" trick, see for instance Section 2.1 
(entitled "The trick") of [5]. It amounts in these cases to the following change of dependent and 
independent variables, 

~ / \ 1 ■ \ / \ exp(io;i) — 1 

z n [t) = ex.p(iojt)z n \T), t = , (2.12a) 

iu 

implying 

z n {0) = z n (0), z n {0) = 4(0) + iuz n (Q). (2.12b) 

Here and below uj is again an arbitrary real constant to which we associate the period T, 
see (2.8), and of course appended primes denote differentiations with respect to the argument 
of the function they are appended to (hence z' n {r) = dz n (T)/dr). Hence clearly the Newtonian 
equations of motion of these two models read as follows: 

z n = 3iojz n + 2u) 2 z n - "SznZn + Siujzl - zi 
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N 

+ E 



2(z n - ibJZ n + CZn)(ze — 10JZ(, + Czj) 



(2.13) 



again with c = 1 respectively c = 1/2. The solutions of the corresponding initial- value problems 
are clearly given by the following variant (of the case with ujj = 0) of Proposition 2.1: 

Proposition 2.2. The solution of the initial-value problems of the two many-body models cha- 
racterized by the Newtonian equations of motion "of goldfish type" (2.13) are given by the eigen- 
values of the N x N matrix U (t) given by the following formulas in terms of the 2N initial 

data z n (0), z n (0) and of the time t: 

U(t) = exp(iut)(l + At + Bt 2 )~ 1 (A + 2Bt), (2.14a) 
with the two constant N x N matrices A and B defined as follows: 

A = diag[5 n (0)], B nm = ^v n v m , v n = z n (0) - iujz n (0) + cz£(0), (2.14b) 

of course always with c = 1 respectively c = 1/2 and r = r(t) defined in terms of the time t 
by (2.12a). Note that the N x N matrix A is diagonal and the N x N matrix B is now dyadic. 

It is plain that the matrix U(t) is isochronous with period T (see (2.14a), (2.12a) and (2.8)), 

U(t±T) = U(t), (2.15) 

and the same property of isochrony holds therefore for the generic solutions of the two many- 
body models (2.13), up to the observation made above (see Remark 2.3). 

3 Proofs 

The starting point of our treatment is the following system of two coupled matrix ODEs satisfied 
by the two N x N matrices U = U(t) and V = V(t): 

U = -U 2 + V, V = -UV + al + 0U + -/V. (3.1) 

Notation 3.1. The 3 scalars a, (3, 7 are 3, a priori arbitrary, constant parameters; I is the 
unit N x N matrix; and we trust the rest of the notation to be self-evident (see also Sections 1 
and 2). 

Remark 3.1. Additional parameters could of course be introduced by scalar shifts or rescalings 
of the dependent variables U, V or of the independent variable t ("time"). We forsake any 
discussion of such trivial transformations (see Remark 2.2). 

To solve this matrix system we introduce the N x N matrix W = W(t) by setting 

U(t) = [W{t)]- l W(t), V(t) = \W(t)\- l W(t). (3.2) 

It is then easily seen that the system (3.1) entails that the matrix W satisfy the following linear 
third-order matrix ordinary differential equation (ODE): 

W = aW + (3W + jW; (3.3) 

and the converse is as well true (in fact, perhaps easier to verify), namely if W satisfies this 
linear third-order matrix ODE, then the two matrices U and V defined by (3.2) satisfy the 
system (3.1). 
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Clearly the general solution of this linear third-order matrix ODE reads 

3 

W{t) = [W ij) exp(iwj-t)] , (3.4) 
i=i 

where the 3 constant matrices W^' are arbitrary, i is the imaginary unit [i? = — 1; introduced 

here for notational convenience), and the 3 scalars ujj are the 3 roots of the following cubic 
equation in u: 

co 3 + i-fu? + /3u - iu = 0, (3.5a) 

w 3 + ijoj 2 + flu — ia = (oo — wi)(u> — w 2 )(a; — CJ3), (3.5b) 

so that 

a = -iojiUJ2UJ3, f} = UJ1UJ2 +^2^3 + 7 = i(u)i + w 2 + w 3 ). (3.5c) 

It is then easily seen that the general solution of the system (3.1) can be written as follows: 

U(t) = i{l + Aexp[i(co2 — ui)t] + B exp[i(u; 3 — u)\)t]} 1 

x \oj\I + w 2 ^4 exp[i(cj 2 — io\)t] + w 3 .B exp[z(w 3 — wi)i]}, (3.6a) 

V(t) = —{I + ^exp[io; 2 - Ux)t] + Sexp[i(w 3 - wi)t]} _1 

x {cjj/ + a;|^4exp[i(cj 2 — + oj$B exp[i(u; 3 — }, (3.6b) 

where A and -B are two, a priori arbitrary, constant N x N matrices. And a trivial if tedious 
computation shows that these formulas provide the solution of the initial-value problem for the 
system (3.1) if the two matrices A and B are expressed in terms of the initial values Uq = U (0), 
Vb = V(0) as follows: 

A = -(wi - w 3 )(w 2 - w 3 ) _1 [l - i(u\ + w 3 )(Vb + wf) 1 (C/b - iw 3 ) _1 ] 

x [j-i( W2 + W3 )(y + w |)" 1 (f/ -ia; 3 )- 1 ]" 1 , (3.7a) 

B = (wi - w 2 )(w 2 - w 3 ) _1 [J - i(wi + <j 2 )(Vb (C7 - «w 2 ) _1 ] 

x [/-i( W2+W3 )(y + W2 2 )-i ([/o _^ 2) -i]-i. (3 . 7b) 

To derive, from the solvable matrix system (3.1), the solvable many-body problem reported 
in the preceding section, we follow the procedure outlined in Section 1. This requires that we 
introduce - in addition to the diagonal N x N matrix Z respectively the nondiagonal N x N 
matrix Y associated to U respectively V via (1.2) respectively (1.3a) - the auxiliary N x N 
matrix M = M(t) defined as follows in terms of the diagonalizing matrix R(t), see (1.2): 

M(t) = [Rit)}' 1 ^). (3.8a) 

In the following we indicate as fi n = [A n (t) respectively M nm = M nm [t) the diagonal respectively 
off-diagonal elements of this matrix: 

M nm = 8 nm /l n + (1 — 5 nm )M nm . (3.8b) 

Remark 3.2. As implied by Remark 1.1, the diagonal elements [i n can be assigned freely, since 
the transformation R(t) R(t) = R(t)D(t) with D(t) = diag[d n (i)] implies, for the diagonal 

elements p, n (t) of the matrix M = [R^)]^ 1 R(t), the expression jl n (t) = fx n {t) + d n (t) / d n (t) , 
with a corresponding change of the off-diagonal elements of the matrix M, M nm =>• M nm = 
(5" 1 M nm 5 m . Note that here and hereafter we denote as M nm the off-diagonal elements of the 
matrix M. 
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It is then easily seen that the equations (1.2) characterizing the time evolution of U and V 
imply the following equations characterizing the time evolution of Z and Y: 

Z + [M,Z} = -Z 2 + Y, Y+[M,Y] = -ZY + aI + f3Z + ~fY. (3.9) 

Notation 3.2. The notation [A, B] denotes the commutator of the two matrices A, B: [A, B] = 
AB - BA. 

Let us now look separately at the diagonal and off-diagonal parts of these two matrix equa- 
tions, (3.9). 

The diagonal part of the first of these two equations reads (see (1.2) and (1.3b)) 

z n = -z n + y n , (3.10a) 
implying 

y n = z n + z 2 n . (3.10b) 
Likewise, the off-diagonal part of the first of these two equations reads 

-(z n - z m )M nm = Y nm , n^m, (3.11a) 
implying 

M nm = n + m. (3.11b) 

Zn Z m 

The diagonal part of the second of these two equations reads (see (1.2) and (1.3b)) 

N 

y n = -z n y n + a + (3z n -\-jy n + ^ (Y ne M en - M ne Y en ), (3.12a) 

implying, via (3.10b) and (3.11b), 

in = -z n z n + jz n - z n + ^z n + pz n + a + 2 ^ ( ^ tn ) ■ (3.12b) 

We now note that, via this equation, time-differentiation of (3.10a) yields the following set 
of Newtonian-like equations of motion: 

z n = -Sz n z n + 7Z n - z n + jz n + /3z n + a + 2 ^ ( YntYin J , (3.13) 

confirming the treatment outlined in Section 1, see in particular (1.4). 

Finally we consider the off- diagonal elements of the second of the matrix equations (3.9). 
The relevant equations read, componentwise, as follows: 

N 

Y nm = -(z n - j)Y nm + y~](Y nk M km ~ M nk Y km ), n ^ m, (3.14a) 

k=l 

namely, via (1.3b) and (3.8b), 

Y n m — \Zn ^Ynm (A*ti \^rn)Y n m ~i~ (.Un Vmi^nm 
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N 

+ Yl (XneMim- M ne Y im ), n^m. (3.14b) 
And via (3.10b) and (3.11b) (and a tiny bit of algebra) this becomes 

Y n m Oil Z " n ^ m 

-ZZ n Z m ~r "f fl n ~r fJ'n 



— it, -ub ' i i — it, ' r^ui 
nm Z n Z m 

N r 

Y n gYt m ( 1 1 



+ E 



Y nm V %l Zn 



n to. (3.14c) 



The next step is to try out the ansdtze (1.5), to see if one can thereby get rid of the quanti- 
ties Y nm . 

We leave the (rather easy but unfortunately unproductive) task to verify that the 
ansatz (1.5a) does not work, i.e. that it does not allow to eliminate the quantities Y nm by 
finding an assignment of the functions g <yl \z) and g^ 2 \z) which, when inserted in (1.5a), yield 
iV(iV — 1) quantities Y nm that satisfy the N(N — 1) ODEs (3.14c) (even by taking into account 
the possibility to assign freely - see Remark 3.2 - the N quantities fi n ). 

We show that instead the ansatz (1.5b) allows the elimination of the quantities Y nm and leads 
to the Newtonian equations of motion "of goldfish type" (2.1). Indeed the insertion in (3.14c) 
of (1.5b) with 

gW{z) = g {2) {z) = g(z), = f®(z) = f(z), (3.15a) 

and the assignment (see Remark 3.2) 

fin = -y (3.15b) 
entails that the N(N — 1) equations (3.14c) can be re- formulated as follows: 



1 \ z n + z 

o \ —■ — , 77 x + —r^ + 3z n ~ 7 + ((n =*-m))> = 

2 [ z n + f(z n ) g(z n ) J z n -z. 



+ Yl {9{zi)[zi + f(zi)](—^ + 1 )}, n^m. (3.16a) 

w.^n.m 1 KZn Z£ Zm ZlJ) 

Notation 3.3. Here and below primes indicate differentiations with respect to the argument of 
the functions they are appended to; and the convenient shorthand notation +((n => m)) denotes 
addition of whatever comes before it, with the index n replaced by the index m. 

It is easily seen that these equations can be re-written as follows: 

Zn + Znf{Zn) + 9 -^4 + 3z n - 7 - Y { 9{zi)[k + m] \ + {{n^m)) (3.16b) 

n/m. 



Zn + f{z n ) g(z n ) iJ^l^n ^ Zn ~ U 

[g(z n ) - l]z n - [g(z m ) - l]z m + g(z n )f(z n ) - g(z m )f(z m ) 

Zn Z m 



This suggests the assignments 

g(z) = l, f(z) = a + bz + cz 2 , (3.17) 
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with the 3 parameters a, b, c a priori arbitrary, since thereby the N(N — 1) equations (3.16b) 
get reduced to the ./V equations 



z n + z n (b + 2cz n ) 



N 



z n + a + bz n + cz\ 
or equivalently (in Newtonian form) 



(2c-3)z„ + 7 + 6 + 2 ^ 



+ a + bz£ + cz| 
z n — Zl 



(3.18a) 



z n = -3z n z n + 7i n + (2c - 3)cz\ + [7c + 3b(c - l)]z n 
+ [o(2c - 3) + 6(7 + 6)]z„ + 0(7 + b) 

2 (i n + a + + cz^) (i^ + a + 6^ + 



N 

+ E 



cz 



Zl 



(3.18b) 



Consistency requires now that this set of N Newtonian equations of motions coincide with 
the N analogous equations (3.13), which, via (1.5b) with (3.15a) and (3.17), now read 



3z n z n -\- ^yz n 

N 

+ E 



2 (i n + a + 5z n + cz^) (i^ + a + bzt + cz|) 



(3.19) 



This clearly requires that the following 4 constraints on the 6 parameters a, j3, 7, a, 6, c 
be satisfied (note that, somewhat miraculously, the two velocity-dependent one-body terms in 
the right-hand sides of the last two equations match automatically, as well as the two-body 
terms): 



(2c-3)c = -1, 
7c + 3b(c - 1) = 7, 
a(2c - 3) + 6(7 + b) = 0, 
0(7 + b) = a. 

And it is easily seen that this entails two alternative possibilities: 
model (i): a, b, 7 arbitrary and 



c=l, j3 = — a + 6(7 + 6), a 
model {it): a, b arbitrary and 

7 = -36, P = -2(a + b 2 ) 



a( 7 + 6); 



1 

2' 



a 



-2ab. 



(3.20a) 
(3.20b) 
(3.20c) 
(3.20d) 



(3.21a) 



(3.21b) 



(Note that, in case (i), another miracle occurred: the solution c = 1 of the first, (3.20a), of the 
4 constraints (3.20) entailed that the second, (3.20b), of these 4 equations hold identically). 

Clearly these two possibilities correspond to the two solvable many-body models "of goldfish 
type" (2.1). 

Next, we must justify the assertions made in the preceding section (see Proposition 2.1) 
concerning the solution of the initial- value problems for the many-body models characterized by 
the Newtonian equations of motion (2.1). The treatment given above (in this section) entails 
that these solutions are provided by the eigenvalues of the N xN matrix U(t) evolving according 
to the explicit formula (3.6a) with (3.7); the missing detail is to express the two initial N x N 
matrices Uq = U(0) and Vq = V(0) appearing in the right-hand side of (3.7) in terms of the 2A^ 
initial data, z n (0) and z n (0), of the many-body problems (2.1). 
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To simplify the derivation of these formulas it is convenient to make the assumption (allowed 
by the treatment given above) that the matrix U(t) be initially diagonal, namely (see (1.2)) that 

R(0) = I. (3.22) 

This entails (see (1.2) and (1.3a)) that 

U(0) = diag[z n (0)], V(0) = Y(0). (3.23) 

The formula (2.3) for Uq is thereby immediately implied. 

The formula (2.4) for Vo = V(0) = Y(0) is also easily obtained, since the expression (3.10b) 
of the diagonal elements of the matrix Y and the ansatz (1.5b) with (3.15a) and (3.17) for the 
off-diagonal elements of Y entail that, componentwise, 

<Wi n (o) + 4(o)] + (i - <W) [z n (0) + a + bz n (0) + cz^(0)] 

x [z m (0) + a + bz m (0) + c4(0)] 1/2 , (3.24) 

and this clearly yields (2.4). 

This completes the proof of the results of the preceding Section 2. As for the findings reported 
in Section 2.1, we consider their derivation sufficiently obvious - for instance via the treatment 
detailed in Section 2.1 of [5]; and by repeating the relevant treatment as given above, especially 
in the last part of this section - to justify us to dispense here from any further elaboration. 

4 Additional findings 

In this section we introduce the time-dependent (monic) polynomial ip(z,t) whose zeros are the 
N eigenvalues z n (t) of the N x N matrix U(t): 

ip(z,t) = deh[zl- U(t)}, (4.1a) 

N N 
${Z, t) = l[[z- Z n (t)} =Z N +J2 [c m (t)z N - m ] . (4.1b) 

?1=1 771 = 1 

The last of these formulas introduces the iV coefficients c m = c m (t) of the monic polyno- 
mial ip(z,t); of course it implies that these coefficients are related to the zeros z n (t) as follows: 

N N 

ci = - y^z n , c 2 = ^2 z n z m , (4-lc) 

n=l n,m=l,n>m 

and so on. 

The fact that the initial-value problem associated with the time evolution (2.1) of the N 
coordinates z n can be solved by algebraic operations implies that the same solvable character 
holds for the time evolution of the monic polynomial ip(z,t) and of the iV coefficients c m (t). 
In this section we display explicitly the equations that characterize these time evolutions. The 
procedure to obtain these equations from the equations of motion (2.1) is quite tedious but 
standard; a key role in this development are the identities reported, for instance, in Appendix A 
of [5] (but there are 2 misprints in these formulas: in equation (A. 8k) the term (N + 1) inside 
the square bracket should instead read (iV — 3); in equation (A. 81) the term TV 2 inside the square 
bracket should instead read A^(A^ — 2)). Here we limit our presentation to reporting the final 
result. 
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The equation characterizing the time evolution of the monic polynomial ip(z,t) implied by 
the Newtonian equations of motion (2.1) reads as follows: 

tpu + {no + m z + V2z 2 )tpzt + (do + &iz)ip t 

+ (ao + ol\Z + c*2Z 2 + a3Z 3 + CX4Z 4 ) ip zz + (/?o + Piz + /3 2 z 2 + f33Z 3 )lpz 

+ {lo + 1iz + j2Z 2 )iJ = 0. (4.2a) 

Here the subscripted variables denote partial differentiations, and the 17 coefficients appearing 
in this equation are defined in terms of the quantities c\ = c\(t), t\ = ci(t) and C2 = C2(t), 
see (4.1c), of the arbitrary positive integer N and of the 3 free parameters a, 6 and 7 charac- 
terizing model (i), respectively of the 2 free parameters a and b characterizing model (ii), as 
follows: 
model (i): 

7] = -2a, 771 = -26, 7] 2 = -2; 

O = -2ci + 2(JV- 1)6-7, 0i = 2JV-l; 

ao = o 2 ) a i = 2a6, «2 = 2a + 6 2 , «3 = 26, 04 = 1; 

/3 = [2ci - (2iV - 3)6 + 7 ]a, A = 26ci - (2JV - 3) (a + 6 2 ) + 67, 

/3 2 = 2ci-2(2iV-3)6 + 7, fa = -(27V - 3); 

70 = -c 1 + c\- 2(N - l)bd + jci -Na + N(N - 2)b 2 - Niry, 

71 = -(2N- l)ci + 2iV(iV _ 2)6 - iV 7 , 7 2 = iV(iV-2); (4.2b) 

model (ii): 

770 = -2a, 771 = -26, 772 = -1; 
O = -ci + 2(N + 1)6, 0i=iV+l; 



«o = o 2 ) ai = 2ab, a 2 = a + b 2 , 03 = 6, 04 



A) = (ci - 2JV6)a, Pi = b Cl - Na-2Nb 2 , (3 2 = ^-2Nb, & = ~; 

70 = -2ci + c 2 - (iV + 2)6ci - ^ + iVa + 7V(iV + 1)6 2 , 

71 = -<™ +J V (JV+1 ) 6 , 72 =^±il. (4.2c, 

Remark 4.1. The equation (4.2a) characterizing the time evolution of the monic polynomial tp 
looks like a linear partial differential equation, but it is in fact a nonlinear functional equation, 
because some of its coefficients depend on the quantities c\ and C2 which themselves depend 
on ifi, indeed clearly (see (4.1)) 

_ ^-P(0,t) _ ^- 2 )(o,t) 

c i = c i(*) = (jV — 1)! ' C2 = C2 ^ = (AT -2)! ' (43a) 

where we used the shorthand notation tp^\z,t) to indicate the j-th partial derivative with 
respect to the variable z of ip(z,t), 

*(»(,,,). **M j = 1 , 2 „... (4 .3 b) 
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As for the system of N nonlinear autonomous second-order ODEs of Newtonian type cha- 
racterizing the time evolution of the ./V coefficients c m , for model (i) they read as follows: 

Cm ~i~ P m ^ C m —1 -|- P m ^ C m ~\~ Prn ^m+l 

+ (- 2 ) c „ + (-l) c 1 + (°) c + flW/j + ( 2 ) C 

= 2cic m - 2(N — m + l)acicm-i 

+ [ci - c 2 + 2(m - l)feci - 7Ci]c m + (2m + l)cic m+ i, (4.4a) 

with the 3JV (time-independent) coefficients pm and the 5iV (time-independent) coefficients qm 
defined here in terms of the 3 free parameters a, 6 and 7 (and of the numbers m and N) as 
follows: 

pfcV = -2(N - m + l)a, p$ = 2(m - 1)6 - 7, p«=2m + l; 
g^; 2 ) = (AT _ m + 2)(JV - m + l)a 2 , 
g*" 1 ) = (JV-m + l)[-(2m-3)6 + 7]o, 

= - m (2N -2m + l)a + m(m - 2)6 2 - m&7, 

= 2(m 2 - 1)6- (m + 1)7, g^=m(m + 2). (4.4b) 
The analogous equations for model (ii) read as follows: 

Cm ~i~ Pm ^C-m—1 Prn Cm "I'Pm^Cm+l 

+ <7m 2) c m -2 + q^Cm-i + 'cm + q$Cm+i + 
= cidm - (JV — m + l)acic m _i 



+ 



2ci + -y - c 2 + (m + 2)6ci 



m 



Cm+ — + 1 ) ClCm+l, (4.5a) 



with the 3JV (time-independent) coefficients pm and the 5iV (time-independent) coefficients qm 
defined here in terms of the 2 free parameters a and 6 (and of the numbers m and N) as follows: 

Pm -V = - 2 (N-m + l)a, p$ = 2(m + 1)6, p« = m + 2; 

= (JV — m + 2)(JV — m + l)a 2 , g*" 1 ) = -2m(N - m + l)a6, 
g^) = —m(N — m — l)a + m(m + 1)6 2 , 

«£> = (m + l)(m + 2)6, ,<? = ( '" + 2) 4 ( "' + 3) . (4.5b) 

Of course in these equations of motion, (4.4a) and (4.5a), it is understood that, for n < 
and for n > N, the coefficients c n vanish identically, c n = 0, while Co = 1 (see (4.1)). 

It is plain that these equations of motion, (4.4) and (4.5), inherit the properties of the original 
many-body models: they clearly are as well solvable by algebraic operations, and of course if 
the original many-body model is isochronous the corresponding model for the coefficients c m is 
as well isochronous, namely 

Cm(t±T) = c m {t), (4.6a) 

and likewise if the original many-body model is asymptotically isochronous, the corresponding 
model for the coefficients c m is as well asymptotically isochronous, namely 

Km[c m (t±T)-c m (t)] = 0. (4.6b) 

t— »oo 
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This of course entails that the conditions on the 3 a priori free parameters a, b and 7 of the 
version (4.4) of model (i), or on the 2 a priori free parameters a and b of the version (4.5) of 
model (n), which are necessary and sufficient to imply that these models be isochronous re- 
spectively asymptotically isochronous, are those indicated in Section 1, namely are those implied 
by (2.6) with (2.5b) and (2.7a) respectively (2.9a); while in the generic case, see (2.10), the 
coefficients c m (t) tend asymptotically to time-independent values, 

lim [c m {t)\ = c m (±oo). (4.6c) 

t— >±oo 

4.1 Two additional isochronous models 

Here we report formulas analogous to those reported above, but related to the many-body models 
discussed in Section 2.1 rather than to those reported in Section 2. We consider the derivation 
of these results sufficiently straightforward not to require any detailed elaboration here beyond 
the terse hints provided below. 

The starting point are the two special cases of the two systems (4.4a) and (4.5a) which obtain 
by setting all the free parameters to vanish. We write them in compact form as follows: 

Cm = 2cc\c m - [2c(m - 1) + 3}c m +i + [(3 - 2c)c\ - c\ + 2(l - c 2 )c 2 ]c m 

+ (2c 2 m + l)cic m+ i - (m + 2) [c 2 {m - 1) + l] c m+2 (4.7) 

with c = 1 respectively c = 1/2. It is remarkable that the number N does not appear in these 
equations; although we always assume that these equations hold for m = 1,2, ... ,N and that 
the dependent variables c n vanish identically for n > N, with N an arbitrary positive integer. 
The diligent reader may also check that this equation also holds identically for m = with cq = 1 
(see (4.1)). 

We then make the following change of dependent and independent variables: 

exp(iwi) — 1 

Cm{t) = exp(imujt)Cm{r), r= — , (4.8) 

with u> again a real arbitrary constant to which we associate the period T (see (2.8)). One thereby 
easily sees that the new dependent variables c m (t) satisfy the following system of autonomous 
Newtonian equations of motion: 

c m = [(2m + l)iu + 2cci]c m - [2c(m - 1) + 3]c m+ i 

+ {m(m + 1)uj 2 - [2c(m - 1) + 3]iwci + (3 - 2c) ci - c\ + 2(l - c 2 )c 2 }c m 
+ {(m + l)[2c(m - 1) + 3]»w + (2c 2 m + l)ci}c m+ i 

-(m + 2)[c 2 (m-l) + l]c m+2 . (4.9) 

It is plain from the way these two models (with c = 1 or c = 1/2) have been derived that 
they are isochronous, namely the generic solutions of these nonlinear Newtonian equations of 
motion satisfy the condition 

Cm(t + T)=C m {t). (4.10) 



5 Outlook 



It is clearly far from trivial that the Newtonian many-body models introduced in this paper - 
see (2.1), (2.13), (4.4), (4.5), and (4.9) - can be solved, for arbitrary initial data, by algebraic 
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operations. Also remarkable is that, for special assignments of their parameters, the systems of 
autonomous nonlinear ODEs (2.1), (4.4) and (4.5) are isochronous or asymptotically isochronous, 
and the systems of autonomous nonlinear ODEs (2.13) and (4.9) are isochronous. 

Let us end this paper by pointing out that Diophantine findings can be obtained from a non- 
linear autonomous isochronous dynamical system by investigating its behavior in the infinitesi- 
mal vicinity of its equilibria. The relevant equations of motion become then generally linear, but 
they of course retain the properties to be autonomous and isochronous. For a system of linear 
autonomous ODEs, the property of isochrony implies that all the eigenvalues of the matrix of 
its coefficients are integer numbers (up to a common rescaling factor). When the linear system 
describes the behavior of a nonlinear autonomous system in the infinitesimal vicinity of its equi- 
libria, these matrices can generally be explicitly computed in terms of the values at equilibrium 
of the dependent variables of the original, nonlinear model. In this manner nontrivial Diophan- 
tine findings and conjectures have been discovered and proposed: see for instance the review 
of such developments in Appendix C (entitled "Diophantine findings and conjectures") of [5]. 
Analogous results obtained by applying this approach to the isochronous systems of autonomous 
nonlinear ODEs introduced above will be reported in a separate paper if they turn out to be 
novel and interesting. 
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